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Abstract 

We suggest a new mechanism for generation of large-scale mag- 
netic field in the hot plasma of early universe which is based on the 
parity violation in weak interactions and depends neither on helicity 
of matter turbulence resulting in the standard a-effect nor on general 
rotation. The mechanism can result in a self-excitation of an (almost) 
uniform cosmological magnetic field. 

The large-scale magnetic field self-excitation in astrophysical bodies like 
Sun, stars, galaxies etc is usually connected with so-called a-effect, i.e. a 
specific term in the Faraday electromotive force £ = aB that is proportional 
to the large-scale magnetic field B. This term is connected with a violation 
of mirror symmetry of a rotating stratified turbulence or convection: the 
number of right-handed vortices systematically differs from the number of 
left-handed vortices due to the Coriolis force action. In this sense, a is 
determined by helicity of turbulent motions. 
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The differential rotation Q usually participates with the a-effect in dy- 
namo action (so-called af^-dynamo), however a-effect induced in a rigidly 
rotating turbulent body could lead to a dynamo action alone (so-called a^- 
dynamo) while the differential rotation alone is unable to result in a dynamo 
action [1]. The a-effect is induced by Coriolis force which destroys the mirror 
symmetry of turbulent motions. 

On the one hand, the a-effect is impossible in electrodynamics of classical 
nonmoving media because of its mirror-symmetry. On the other hand, the 
mirror asymmetry of the matter happens at the level of particle physics and 
we can expect that an a-cffect could be based on this asymmetry. The aim 
of our paper is to present such mechanism based on parity violation in weak 
interactions. 

We remind that the main problem of most particle physics mechanisms 
of the origin of seed fields is how to produce them coherently on cosmolog- 
ical (large) scales. There are many ways allowing to generate seed small- 
scale random magnetic fields in early universe, e.g. at phase transitions [2], 
however, the following growth of the correlation length, e.g. in the inverse 
cascade with the merging of such small-scale fields [3] , hardly could produce 
a substantial large-scale fields at present time [4]. We do not consider an 
evolution of correlated domains and the corresponding growth of correlation 
length considered e.g. in the review [5] and concentrate here on the gener- 
ation of a mean magnetic field (amplification of its strength) via a-effect if 
such mean field has been already seeded somehow from small-scale magnetic 
fields. 

Let us consider hot plasma of early universe after electroweak phase tran- 
sition, T ^ Tew — 10^ MeV, when we may use point-like (Fermi) approx- 
imation for weak interactions and where at the beginning a weak random 
magnetic field has a small macroscopic scale comparing with the horizon, 
A<^ Ih, while within a domain of the volume ~ such magnetic field can 
be uniform and directed along an arbitrary z-axis, B = (0, 0, B). Obviously, 
this docs not violate the isotropy of universe as a whole with many randomly 
oriented domains. 

Within a domain with an uniform magnetic field obeying the WKB limit 
I e I S <S the single quantum (spin) effect remains for electrons and 
positrons which populate the main Landau level only and contribute to the 
lepton gas magnetization, Mj'^^ = fis < i'aljlbi'o- >= ^jzl^B{^&^ o")^Ocr ~ 
(sgncr)i? [6], where [ib is the Bohr magneton, noo- is the number density at 
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the main Landau level for the electrons and positrons (cr = e , cr = e"*"), 

\e \ B , , \e\B TIn2 

The magnetization Mj^-* changes sign for electrons and positrons , (sgn cr) = 

±1, effectively due to the opposite spin projections on the magnetic field at 
the main Landau levels. 

For a small magnitude of magnetic fields we neglect small polarization 
of other components: muons, tau-leptons, quarks or nucleons. Obviously, 
densities (1) are very small comparing the total lepton densities 

= /(rfV(2vr)=')/i?(%) ~ 0.183 rio. < n,. Here a = e^; is 
the Fermi distribution; e, £p = \/p^~-i-m|, T are the lepton electric charge, 
the energy and the temperature of lepton gas correspondingly. 

In magnetized plasma the pseudovector Mj'^^ enters weak interaction of 
the charged a-fiuid component with neutrinos (antincutrinos) through the 
axial part of the point-like current x current interaction Hamiltonian, 



where Gp — 10"^ /rrip is the Fermi constant, ^j^'^^ — ju-ju is the neutrino 
current density asymmetry. 

Such interaction provides a force F"'''*'^ (see below Eq. (6)) that is addi- 
tive to the Lorentz force qa(Ei + V^. x B) acting in MHD plasma on charged 
particles of the kind a and obviously depends on gradients of the interac- 
tion potential, (F^^^)i ~ —diVj'^\ or for an uniform magnetization within 
a domain (Mj = constant) on derivatives of the neutrino current density 
asymmetry, di5j^p . 

The electric field E being common for all charged particles is obtained 
multiplying the motion equations for each charged components by the cor- 
responding electric charge with the following summing over components 
that leads to Yl,a Lorentz force and the remarkable addition of 

noo- in the weak electromotive force term for electron-positron components, 
j^jcak _ ^j^, ^ {uq^ + ^0-), due to the independence of the product 

qadiVa"^^ ~ qfjGpMj'^^ ■ didj'^p / jis on the sign of the electric charge since 
Mj-'^'* ~ (sgn a) . Such term violates parity and provides a new particle physics 
origin of a-effect for magnetic field generation, d^Q = —V x E™^*^^. 
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The pair motion equation in the one-component MHD is derived after 
the summation of Euler equations for comoving electrons and positrons for 
which the standard (polar vector) electric field cancels since — ± \ e \, the 
standard Lorentz force | e | (V+ — V_) x B = rotB x B/47rne arises while 
the weak force term depends in hot lepton plasma on the negligible difference 
of densities, (no+ — ^o-) as well as the neutrino axial vector potential V^^^ 
describing a probe neutrino in the electron-positron plasma [7] when 6jj^^ — >■ 
kj and the sum over a leads to V^^) = G'i.^(no+ - no-)B ■ k/Bk. 

Now we estimate the a-effect originated in early universe by particle 
physics effects. Let us note that in an external large-scale magnetic field 
B a polarized equilibrium lepton plasma is characterized by the density ma- 
trix, 

/^^^ (e,) = %4^) is,) + ^^e? is,) , (2) 

where a is the Pauli matrix; b = B/i? is the ort directed along the magnetic 
field; feq{£p) is the Fermi distribution; siq{ep) = — (| e | B /2ep)dfeq \£p) /dsp 
is the spin equilibrium distribution that defines the number density at the 
main Landau level (1), / d^pS^\ep) / (27r)^ — noa', A = ±1 is the spin pro- 
jection on magnetic field. 

Then we start from the linearized relativistic kinetic equations (RKE) 
derived in Vlasov approximation for a magnetized lepton plasma in the stan- 
dard model (SM) of electroweak interactions after the summing over spin 
variables as given in Eq. (30) of Ref . [6] , 

dt dx 

op op 

'^"-^weak "'"'weak ^p ~ \'^) 

Here weak forces F^^j^, F^^j^ that appear due to the generalization in SM 
of the standard Boltzman equation have the form: 



f2:L= (sgna)GpV^^civ) 
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+v X V X 5j^J , 



(4) 



and 



Fl1L= -(sgna)GpV^^c(^) 



-V X V X 5n^^h + — V(a(p) • Sj^J 



dt 

(5) 



c^^ — 2^±0.5, ci'^'* = =f0.5 are the vector and axial couplings correspondingly 
(upper sign for electron neutrino) where subindex a = e, fi,T characterizes 
the kind of neutrino, ^ = sin^^vi^ — 0.23 is the Weinberg parameter; Sj^^ — 
j^a ~ -^^a neutrino four-current density asymmetry, 

/d^k - -J 

(27r)^ " ^) 

is the neutrino (antineutrino) four-current density; 5ni,^ — n^^ — is the 
neutrino density asymmetry that plays an important role in the generation 
of magnetic field (see below (10) and in [8]). 

Finally a(p) in the last term in (5) is the three-vector component of the 
four- vector that is the analogue of the Pauli-Lubahski four- vector a^ip) = 

Tr (p757^)/4me = [pC/me] ( + p{p ■ C)/me{Sp + meij with the change of 

the spin ( to b. 

Notice that we can substitute the total number density distribution func- 
tion /^'^^(p, X, t) — feq\sp) + Sf^'^\p, X, t) normalized on the total density 
n„ = f {d^p/{27rf)f^''\p, X, t) / (rfV(27r)^)/i? into all the terms of RKE 
in the first and second lines (3) restoring its standard Boltzman form. 

In addition, RKE (3) obeys the continuity equation, dj\^ ^ /dx^ = 0, or the 
lepton number is conserved, where t) = J d^p{pfj_/ Sp) f^'^\p, x, t)/(27r)^ 

is the lepton four-current density. 

Then we can use the standard method [9] for transition from kinetic equa- 
tions to the hydrodynamical ones. Multiplying RKE (3) by the momentum 
p and integrating it over d^p with the use of the standard definitions of the 
fluid velocity Vo- = J d^pvf'^'^\p. x, t) / {2ivY and the generalized momen- 
tum Pp. = J (i^pp/'^'^^(p, X, t)/(27r)^ one obtains the Euler equation for 
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the fluid species a in plasma with the additive collision terms taken in the 
r-approximation, 



a emrn / i \T» ^Pcr , 
-T— = Oi'a - [V^v + Vavjt'a h 

dt ria 

+g,(EH-[V,xB]) + F--'^ . (6) 

Here a = e~, t~ , q^, Qd, ■■■ (o' = e"*", /i"*", r"*", g„, qa, ... for antiparticlcs) enu- 
merates the plasma components; z/""^ is the electromagnetic collision fre- 
quency which leads to the fast equilibrium in plasma and defines plasma 
conductivity; u^^, u^p are the weak collision frequencies providing the gen- 
eration mechanism suggested in [8] for neutrino scattering off electrons and 
positrons before neutrino decoupling; p„ is the fractional pressure. 

We isolate in the weak ponderomotive force vector and axial parts, i.e. 
pwcak ^ pCV') _|_ pC^)^ 'pj^g gj^g^ term Fi^-* coming from (4) was found by the 

independent (Lagrangian) method in [10] and is irrelevant to the magnetic 
field generation mechanism considered here. The axial vector force f'^J^^ 
appearing from (5) due to the polarization of lepton gas in a weak magnetic 
field B, is 



dt 



a=e,iJ,,T 

+Aro,V(b-5jj]. (7) 
Finally the relativistic polarization term A^oo-, 

9(2vr)3 I fk\'r)dp-[v{3-v)], (8) 

in the non-rclativistic case tends to the lepton density at the main Landau 
level given by Eq. (1), Nq^ tiqct- Obviously, the weak force (7) changes 
sign for positrons, cr — > a, due to the signature function. 

For the hot plasma multiplying the Eulcr equation (6) by the electric 
charge q„. summing over a and dividing by XIo-^ct ~ ^"^^ electric 

field E including all known polar vector terms phis the new axial vector 
E ~ aB originated by electron-positron polarizations which violates parity. 
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This is similar to the derivation of E in [10] for unpolarized plasma, and , in 
particular, from the Lorentz force one obtains the term — 'Yliai^^l Q"^)^ cr x B 
that obviously leads from the Maxwell equation 9tB = —V x E to the dynamo 
effect in Faradey equation. Thus, we arrive to a governing equation for 
magnetic field evolution 



where wc omitted dynamo term neglecting any macroscopic rotation in plasma 
of early universe ^ . 

Here we approximate the tensor aij coming in E from the axial vector 
force (7) by the first diagonal (~ a5ij) term: 



where densities no± are given by Eq. (1), equilibrium densities obey n^/rie — 
0.5, and we assume a scale of neutrino fiuid inhomogeneity t ~ Ag^^^j, that is 
small comparing with a large A-scale of the mean mag nctic field B, X^^^.^^ < A. 
Let us stress that the addition of positron and electron contributions in a 
stems from the change of the sign in the weak force (7). 

The origin of the sca/ar a-coefficient (10) from weak interactions leads to 
its important difference from the pseudoscalar coefUcient a —< v-(V x v) > in 
the standard MHD based on fluid electrodynamics, where C, P, T symmetries 
are conserved separately. 

Really, while in the last case (standard MHD) all terms in induction 
equation are pseudovectors obeying PBP~^ = B, etc., in our case the flrst 
term in the r.h.s of Eq. (9) turns out to be a pure vector, violating parity, 
P(V X aB)p-^ = -V X aB. 

Nevertheless, all terms in the generalized induction equation (9) obey CP- 
invariancc as it should be for clectroweak interactions in SM since the new 

^The detailed derivation of full Faradey equation including dynamo, Bierman battery 
terms, both vector weak interaction cotribution [10] and new axial vector interaction terms 
is given also in [11] where lepton MHD equations in SM for unpolarized and polarized 
media are derived from RKE's [6]. 



— = V X aB + T^V^B , 



(9) 




(10) 
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coefficient a (10) is CP-odd, CPa{CP) ^ = —a, as well as curl-operator 
(V X ...). This is due to the changes Uq Ug^'' and ^n^^ — >• —5ny^ in 
(10), provided by the well-known properties: particle helicities are P-odd 
and particles become antiparticles under charge conjugation operation C, 
in particular, the active left-handed neutrinos convert into the active right- 
handed antineutrinos under CP- operation, z/^ — > i>a. 

Finally the diffusion coefficient 7] ~ (47r x 137 T)~^ is given by the rela- 
tivistic plasma conductivity. We do not present in Eq. 9 standard terms like 
differential rotation etc which seems to be not very important in the problem 
under consideration. 

This is our main result. We stress that the Eq. 9 is the usual equation for 
mean magnetic field evolution (see e.g. [12]) with a-effect based on particle 
effects rather on the averaging of turbulent pulsations. It is well-known (see 
e.g. [1]) that Eq. 9 describes a self-excitation of a magnetic field with the 
spatial scale A ~ ry/o; and the growth rate a'^/Arj. 

Let us estimate these values for the early universe. For a small neutrino 
chemical potential /ij,, Cua{T) = iJ,ua{T)/T <^ 1, the neutrino asymmetry in 
the r.h.s. of Eq. (10) is the algebraic sum following the sign of the axial 
couphng, dtl = ±0.5, 

^-E«^-|;)K.m+6.(T)-6.(T)i. (u) 

Wc take for crude estimations below C,Uf, (T) + (T) — (T) ~ ~'^^u^ {T) 
because different chemical potentials almost compensate each other for high 
temperatures [13], i.e. {,,(T) + ^.^(T) + ^..(T) ^ 0. 

As a result, we arrive to the following estimate of the a-coefficient (10), 

a = 2.8 x 10-3^(r/MeV)6(l,(T)/A(;),) | I, 

where a free parameter for our coUisionless mechanism - scale Ag^^j^j is nor- 
mahzed on the neutrino mean free path li,{T) — given by the weak rate 
Tw = 5.54 X 10-22(T/MeV)5 MeV. 

Substituting a into A — r]/ a we arrive now to the estimate 
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where Ih{T) = {2H)-^ and H = 4.46 x 10~'^^{T/MeVf MeV is the Hubble 
parameter. 

If the neutrino fluid inhomogeneity scale Xq^^^ is of the order lt,{To) ~ 
4 cm -C Ih{To) ~ 10^ cm, we have I^/Ir > 1 at the beginning of the lepton 
era (T = Tq ~ 10^ MeV, redshift 2; ~ 3 x 10^^), or more correctly, accounting 
for the BBN hmit ^^^(T) < 0.07 [13] obtained for Tbbn = 0.1 MeV , the 
mean magnetic field will be uniform in whole universe, A//j^ > 1, at T ~ 118 
MeV. If this neutrino parameter would be much smaller at high temperatures 
To, ^;/e(7o) <^ 0.07, one can choose another free neutrino parameter Ag^^^j <S 
li,{To) in such a way that the ratio Agy^;^/(/,y(To) | Cu^iTo) |) remains invariant 
and our conclusion about the tendency to a global uniform field is still valid. 
Note that for the neutrino gas the macroscopic parameter A^^^j^ varies in a 



where -Bmax is some seed value at the instant T^ax Tew ~ lOOGeV (here 
we imbed the standard estimates of a^-dynamo into the context of expanding 
Universe) . 

For Ag'J^^jj(T) ~ LiT) we can estimate the index in the exponent (13) sub- 
stituting in the integrand the expansion time 

t{T)= 3.84xl02^(T/MeV)"2MeV~Vv/g^ with the effective number of degrees 
of freedom g* ~ 100 at the temperatures T > 1 GeV. Then from our esti- 
mates of a{T) and r]{T) with the change of the variable (T/2 • lO^'MeV) — > x 
one finds the fast growth of the mean field (13) in hot plasma, x < 1 with a 
conservative estimate. 



given by the upper limit Xmax = Ij Tmax = 20GeV. Such upper limit defines 
entirely the magnetic field amplification due to the steep dependence on 
the temperature and still obeys the point-like Fermi interaction we rely on. 
As in the case of magnetic field scales (12) the second free parameter Ag^^j^ 
can be chosen much smaller, Xq^^^ <^ ^i^(^) ) providing the invariant ratio 



wide region T'^ < X^^d < U^)- 

The magnetic field time evolution is given by 




(13) 




(14) 
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lu{T) I I /^fiuid ^'^^ very small neutrino chemical potential ^ueiT) ^ 0.07 
and resulting in an enhancement of a small mean magnetic field -Bmax 
^max/ 1^1*^ ^Iw/ I 6 I by collective neutrino-plasma interactions considered 
here in Eq. (14). 

Note that the inflation mechanism (with a charged scalar fleld fluctuations 
at super-horizon scales) explains the origin of mean fleld at cosmological 
scales, however, the value of this fleld is too small for seeding the galactic 
magnetic fields [14]. 

The amplification mechanism suggested in our paper can improve this 
very low estimate by a substantial factor from Eq. (14). 

Thus, while in the temperature region Tew ^ T ^ Tq = 10^ MeV 
there are many small random magnetic field domains, a mean magnetic field 
turns out to be developed into the uniform global magnetic field. The global 
magnetic field can be small enough to preserve the observed isotropy of cos- 
mological model [15] while strong enough to be interesting as a seed for 
galactic magnetic fields. This scenario was intensively discussed by experts 
in galactic magnetism [16], however until now no viable origin for the global 
magnetic field has been suggested. We believe that the a^-dynamo based on 
the CK-effect induced by particle physics solves this fundamental problem and 
opens a new and important option in galactic magnetism. 
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